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Newsletter  No.  10 


Category  III  Trigonometry  Modules 


February  1,  1971 


We  are  enclosing  four  trigonometry  modules.  The  module  "CIRCULAR 
FUNCTIONS"  was  written  by  Mr.  Robert  Gillespie  of  Allderdice.  Its  place 
in  mathematics  is  described  in  the  overall  plan  he  and  his  group  prepared 
last  summer.  The  "12th  grade  mathematics"  topics  given  in  Newsletter  No. 

7 came  out  of  that  plan. 

4 

As  a result  of  mathematical  questions  raised  by  "CIRCULAR  FUNC- 
TIONS", two  other  modules  were  prepared  called  "TRIGONOMETRY  FUNC- 
TIONS AND  TCHEBYCHEV  POLYNOMIALS"  and  "INVERSE  CIRCULAR 
FUNCTIONS".  The  module  "VORTAC"  is  for  the  more  ambitious'  student. 
Students  who  want  to  try  a two  aircraft  simulation  of  navigating  by  VORTAC 
should  type: 

>RUN  166FW  /WGAME/ 


Call  For  Student  Programs 

We  would  like  to  circulate  student  programs  of  all  levels  through  the 
newsletter.  This  is  not  a request  for  "best  programs",  but  rather  a mech- 
anism through  which  teachers  can  give  recognition  to  experienced  as  well 
as  neophyte  student-users  of  computing.  Please  include  a listing,  sample 
run,  and  appropriate  documentation,  including  teacher's  name,  subject, 
student's  name,  age,  and  grade. 
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Exchange  Programs 

We  are  looking  into  the  possibility  of  exchanging  newsletters  and  pro- 
grams with  other  high  sc  ho  o’  computing  efforts.  The  two  newsletters  we  are 
aware  of  right  now  come  from  the  Southern  Minnesota  School  District  Project 
at  Mankato,  Minneapolis,  and  from  Brooklyn  Polytech  in  New  York  (The 
Huntington  Project).  This  later  group  is  specializing  in  simulation  (type  II) 
programs. 


^Supported  in  part  by  NSF  grant  GJ  1077 
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CIRCULAR  FUNCTIONS 


Consider  an  angle  A with  its  initial  side  placed  along  the 
X-axis,  with  its  vertex  at  the  origin  P:(0,0),  and  with  a terminal 
side  that  passes  through  the  point  P:(X,Y).  Let  us  call  the  dis- 
tance from  (0,0)  to  (X,Y)  R,  so  that  R^  = x^  + Y^*  Then  the  six 

circular  functions  of  the  angle  A are  defined  as: 


Sine  of  the  angle  = y/r 
Cosine  of  the  angle  - x/r 
Tangent  of  the  angle  = y/x 


Con tangent  of  the  angle  = x/y 
Secant  of  the  angle  = r/x 
Cosecant  of  the  angle  ~ r/y 


Sample  problem: 


Find  the  circular  functions  of  an  angle  (A)  whose  terminal 
side  passes  through  the  point  P:  (5,12). 


Solution: 

x = 

5 and  y 

therefore 

r ss 

SQRT  (5*5 

thus 

sin 

cos 

A = 12/13 

A - 5/13 

tan 

A = 12/5 

cot 

A = 5/12 

sec 

CSC 

A = 13/5 

A = 13/12 

* 12 

+ 12*12)  * SQRT (169) 

- .9231 

- .3846 

- 2.400 

* .4167 

= -2,600 
= 1.0833 


Using  the  Computer 


On  the  next  page  we  will  show  the  flow  chart  for  a program 
that  asks  the  user  to  supply  values  for  X and  Y.  The  program 
then  calculates  SIN  A,  COS  A,  and  TAN  A. 

Notice  that  care  is  taken  to  avoid  dividing  by  zero.  (For 
which  circular  functions  can  this  happen???) 

More  advanced  programming  problems  are  suggested  on  pages 
5 and  6. 


♦NOTE:  When  using  the  computer  you  should  type  5 as  5.  or  5.0 
and  12  as  12*  or  12.0-  etc. 
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Flowchart  for  a Program  to  Calculate  Sin,  Cos,  and  Tan 


4 


i 


1 


CIRCULAR  FUNCTIONS 


4 


* 


I 


Here  is  an  NBS  program  based  on  the  flow  chart  on  page  3: 

50  PR.  "THIS  PROGRAM  CALCULATES  THE  SIN,  COS,  AND  TAN  FUNCTIONS" 

60  PR.  "OP  AN  ANGLE  A DEFINED  BY  THE  RATIOS  Y/R,  X/R,  AND  Y/X  " 

70  PR.  "RESPECTIVELY" 

95  PR. 

100  PR.  "TYPE  IN  THE  COORDINATES  OF  YOUR  POINT  IN  X,Y  ORDER,  USING" 
110  PR.  "DECIMAL  FORM" 

120  INPUT  X,Y 
200  RfSQRT  (X*X+Y*Y) 

300  S*Y/R 
400  C=X/R 

450  IF  X=0  GOTO  900 
500  T=Y/X 

600  PR.  "SIN(A)=":S?"COS(A)=":C;"TAN(A)=":T 
700  PR.  "DO  YOU  WISH  TO  ENTER  ANOTHER  POINT"? 

710  INPUT  R$ 

720  Y$="YES, SURE, OK, O.K., CERTAINLY, YUP,OF  COURSE, RIGHT, YEAH, Y, " 

730  IF  IEQIV (R$ , Y$ , 1)  GOTO  100 
800  STOP 

900  PR.  "SIN  (A)  : S ? "COS  (A)  ="  s C ; "TAN  UNDEFINED" 

910  GOTO  700 
>RUN 

THIS  PROGRAM  CALCULATES  THE  SIN,  COS,  AND  TAN  FUNCTIONS 
OF  AN  ANGLE  A DEFINED  BY  THE  RATIOS  Y/R,  X/R,  AND  Y/X 
RESPECTIVELY. 

TYPE  IN  THE  COORDINATES  OF  YOUR  POINT  IN  X,Y  ORDER,  USING 
DECIMAL  FORM. 


d?5.0,  12.0 


SIN (A)— .9231  COS (A)=. 3846  TAN (A) =2. 4 

DO  YOU  WISH  TO  ENTER  ANOTHER  POINT?  SURE 
TYPE  IN ETC. 


‘"NOTE:  You  must  type  in  two  decimal  numbers-  If  you  don’t  know  a 
number  (say  the  square  root  of  3) , you  can  use  @NBS. 

EXAMPLE:  Suppose  you  want  to  find  the  circular  functions  for 
P:(/3f  , /2  ). 

>RUN 


WRONG!! 


TYPE  IN  THE  COORDINATES  OF  YOUR  POINT 
*?SQRT(3),  SQRT  (2) 


>RUN 

THIS  PROGRAM  ETC. 


CORRECT 


TYPE  IN  THE  COORDINATES  OF  YOUR  POINT  IN  X,Y  ORDER, 
?@NBS 

>PR.  SQRT (3) 

1.732 

>PR.  SQRT (2) 

1.414 


>EXIT 

4 71.732,  1.414 
*' ETC 


0 


1 
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Computer'  problems 


1.  Use  the  program  on  Page  4 for  P:(l,l),  P:(0,1),  P:(1,0), 

P:(-3,3),  P:(-2,3),  P:  (-3,-3) , P:(3,-3)  £IMPT:  Type  1 as  1.0  etc.] 

2.  Modify  the  program  given  on  Page  4 so  that  all  six  circular 
functions  are  calculated. 

a.  Use  your  computer  program  to  find  the  circular  functions  of 
the  following: 

i.  An  angle  A whose  terminal  side  passes  through 
P: (2*SQRT(2) , SQRT(7>)  NOTE:  You  must  use  @NBS 
to  find  2*SQRT(2)  etc. 

ii.  An  angle  A whose  terminal  side  passes  through 
P:  (-2.386,  7.590) 

iii.  An  angle  A whose  terminal  side  passes  through  P: (0,2.7) 

iv.  An  angle  A whose  terminal  side  passes  through  Ps (2.7,0) 

b.  Run  your  program  again , supplying  data  on  several  points  in 
all  four  quadrants.  Can  you  make  a general  rule  about  the 
signs  of  the  circular  functions  in  each  quadrant? 

c.  Investigate  angles  whose  terminal  sides  are  coincident  with 
the  coordinate  axes; that  is, circular  functions  of  0,  90, 

180,  360°.  Which  functions  are  undefined? 


Advanced  Problems: 

3.  Write  a program  that  automatically  generates  the  six  circular 
functions 

FOR  X — 1.0,  .9,  .8,  .7,  .6,  .5,  .4,  .3,  .2,  .1,  0. 

— .1,  -.2,  -.3,  “.4,  -.5,  —.6,  “.7,  -.8,  ™«9,  —1.0 

WITH  R FIXED  AT  R = 1 (UNIT  CIRCLE)  . 

Here  is  part  of  a program  to  do  this  with  some  output: 

100  AS=*'2<SD.DDB)  2CSD.DDDB)  3< SDD.DDDB)  SDD.DDDB/'* 

no  B$s'*2(SD.DDB)  2<SD. DDDB)  * UNDEF  * SDD.DDDB  * UNDEF  ' SDD.DDDB/*' 
120  CS="2CSD.DDB)  2CSD.DDDB)  SDD.DDDB  ' UNDEF  * SDD.DDDB*  UNDEF  */** 
130  R*1 

140  PRINT  '*  X Y SIN  COS  TAN  COT  SEC  CSC'* 

150  FOR  X=1  BY  -. 1 WHILE  X>-1.1 
160  IF  ABS<X)>M  GOTO  280 
170  Y=ABS<SQRT<R*R-X*X>> 

1 8.0  IF  ABSCXX.005  GOTO  250 
190  IF  ABS<YX.005  GOTO  280 

200  PRINT  IN  FORM  AS*X* Y* Y/R:X/R: Y/X*X/Y*R/X*R/Y 
210  NEXT  X 
230  END 

O 

o 


CIRCULAR  FUNCTIONS 


6 


3.  (Continued) 

>RUN 

X 

Y 

SIN 

C0S 

TAN 

COT 

SEC 

CSC 

+ 1*00 

+0*00 

+0*000 

+1 *000 

+00*000 

UNDEF 

+01*000 

UNDEF 

+ 0*90 

+0*44 

+0*436 

+0*900 

+00*484 

+02*065 

+01*111 

+02*294 

+ 0*80 

+ 0.60 

+0*600 

+0*800 

+00*750 

+01*333 

+01*250 

+01*667 

+ 0*70 

+0*71 

+0*714 

+0*700 

+01*020 

+00*980 

+01*429 

+01*400 

+ 0*60 

+ 0*80 

+0*800 

+0*600 

+01.333 

+00*750 

+01*667 

+01 *250 

+ 0*50 

+ 0*87' 

+0*866 

+0*500 

+01.732 

+00*577 

+02*000 

+01 .155 

+ 0*40 

+0*92 

+0*917 

+0*400 

+02.29J 

+00*436 

+02*500 

+01*09 1 

+ 0*30 

+ 0*95 

+0*954 

+0*300 

+03*180 

+00*314 

+03*333 

+01 

+ 0*20 

+ 0*98 

+0*980 

+0*200 

+04*899 

+00*204 

+05*000 

+01*021 

+ 0*  10 

+0*99 

+0*995 

+0* 100 

+09*950 

+00* 101 

+10*000 

+01*005 

+ 0*00 

+ 1*00 

+1*000 

+0*000 

UNDEF 

+00*000 

UNDEF 

+01*000 

-0*  10 

+ 0*99 

+0*995 

-0* 100 

-09*950 

-00*101 

-10*000 

+01*005 

ETC* 


4 . Mem  Against  Machine: 

Try  to  make  the  computer  blow  its  mind  by  calculating 
TAN (A),  with  R * 1,  and  X approaching  zero* 

10  LET  X=1 

20  IF  X<*0000000000001  G0T0  60 
ONE  30  print  "X="*X*TAB<35>fTAN<X>="sSGRTO-X*X>/X 

APPROACH  40  let  X=X/2 
50  G0T0  20 
60  END 


5*  The  picture  on  the  cover  shows  that  the  paint  dripping  from 
the  paint  dispensers  on  the  three  wheels  (rotated  by  monkey 
power)  traces  out  "graphs"  of  the  circular  functions  on  the 
canvas  (moved  by  mad  artist  power) . 


Write  a program  that: 


(a)  Graphs  the  SIN  or  COS  function* 

OR  (b)  Graphs  the  SIN  and  COS  function  on  the  same  axis,  using 
two  different  values  of  R* 

OR  (c)  Graphs  SIN  and  COS  functions  with  different  values  of  R 
and  different  starting  points,  possibly  making  "computer 
generated"  art  as  your  goal* 

* * 4 

* 1 
Suggestions : 

1.  Ask  your  teacher  for  the  module  called  "Graphing  by  Computer* 

Could' you  make  your  graph  print  symbols  like  "2",  "3",  etc* 
to  make  a "number"  painting  for  later  water  coloring,  where 
1 = Red,  2 = Blue,  etc* 


2* 


inverse  circular 


But  the  principal  value  of  ARCSIN(0,8)  is 
0.9273  RADIANS  WHICH  CAN  BE  CALCULATED  BY 

• ARCS  I N (X)  ^ TT/2  - U 1 - X_ ) * 
(1,570795207-  ,214512362  * X 
+ ,087876311  * X2-  ,044958884  * X* 

+ .019349939  * X*-  ,004337769s*  Xs) 


PROJECT  SOLO 
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FINDING  THE  INVERSE  CIRCULAR  FUNCTIONS 


In  the  module  on  trigonometric  functions  we  saw  that 
polynomial  approximations  cure  used  by  computers  to  calculate 
the  circular  (or  trigonometric)  functions  when  the  indepen- 
dent variable  is  the  angle. The  simple  NBS  program: 


>10  LET  A=.9273 

>20  PRINT  "A=":A;"SIN(A)=":SIN(A) 

>RUN 

will  result  in  the  output: 

A=.9273  SIN(A)=. 800002869 

In  order  that  the  computer  successfully  execute  this  program 
it  must  use  one  of  the  Hastings  approximations  for  finding 
SIN(.9273).  (See  Problem  3 in  the  module  "TRIG.  FUNCTIONS 
AND  TCHEBYCHEV  APPROXIMATIONS".) 

THE  INVERSE  PROBLEM 

Now  suppose  we  ask  the  computer  to 

>PR.  ARCS  IN  (.8) 

we  will  obtain: 

0.927295218 

We  have  asked  the  computer  to  print: 

"the  angle  which  has.  a sine  of  0.8" 
or  to  be  exact: 

"Print  the  angle  A (in  radians)  such 
that  SIN  (A)—  0.8,  and  A lies  in  the 
interval  -ir/2<A£ir/2 

The  cover  of  this  module  shows  why  we  have  to  agree  on  an 
interval  for  A.  The  inverse  function  ARCS IN  is  not  s inale 
valued , and  so  we  must  eliminate  confusion  by  using  the 
(unspoken)  assumption  that  the  principal  value  of  A will  be 
calculated.  For  the  SINE  function,  the  principal  value  is 
found  in  either  quadrant  I or  quadrant  IV. 


HASTINGS  SHEET  35 


The  formula  for  approximating  ARCS IN (X)  shown  on  the 
front  cover  is  Hastings'  Sheet  37,  a rather  accurate  for- 
mula* It  assumes  that  X is  limited  to  the  range 

0<X<1 

Here  is  an  NBS  program  to  use  a simpler  but  less  ac- 
curate formula  (Hastings'  Sheet  35)  to  calculate  the  ARCSIN 
function*  The  program  also  prints  the  more  accurate  NBS 
ARCSIN  function  for  comparison*  The  Sheet  35  formula  is: 

ARCSIN  (X)^ir/2-  CITOC) * (1*  5707288- * 2X21144*X+*  0742610*X*X 
-♦0187293*X*X*X) 

Notice  that  polynomial  is  evaluated  using  nested  multiplica 
tion  in  line  560  in  order  to  reduce  the  nunlber  of  multipli- 
cations* Why  is  this  important? 


PROGRAM  ,0  COMPARE  NBS  LIBRARY  ARCSIN  FUNCTION 
WITH  HASTINGS  35  FUNCTION 


ARGUMENT  NBS  LIBRARY  HASTINGS  35 


0 

0*  1 
0*2 
0*3 
0*4 
0*5 
0.6 
0*7 
0*8 
0*9 
1 

>LI STNH 

500  PREPROGRAM  TO  COMPARE  NBS  LIBRARY  ARCSIN  FUNCTION  " 

501  PREWITH  HASTINGS  *35  FUNCTION*'  PR*  PR* 

510  LET  A0=1  *5707288  LET  Al=-*2121144 

520  LET  A2«* 07426 10  LET  A3*-*0187293 
530  P=3* 14159/2 

5 40  PR E ARGUMENT" :TAB< 20) : "NBS  L I BRARY” : TAB  < 40  it  "HASTINGS  35" 

545  PR." " 

550  FOR  X«0  TO  1 STEP  0*1 
5 60  D=«A3*X+A2)*X+A1)*X  + A0 
570  S=P-S©RTC!-X>*D 

580  PR*X: TAB(20) t ARCSINCX ) : TAB( 40) : S 
590  NEXT  X 
600  END 


0*100167421 
0*201357921 
0*304692654 
0*41 1516846 
0*523598776 
0*643501209 
0*775397497 
0*927295218 
1 * 1 19769515 
1 *570791655 


6*  620000204E-05 
0*  100107027 
0*201313803 
0*304700575 
0.41 1559864 
0*523643617 
0*643521213 
0*775385477 
n 0*927265853 
1*119751877 
1*570789605 


1 


4 


Pr  ob  lean  1 

Improve  the  above  program  so  that  X may  lie  in  the 
interval  -1<X<1.  Test  your  program  by  running  it  with  line 
550  changed  to: 

550  FOR  X=-l  TO  1 STEP  0.1 
[HINT : RECALL  THAT  SIN (-X) =“SIN (X) ] 

Problem  2 

Find  out  what  NBS  does  when  you  type 
>PR.  ARCSIN (2) 

Can  you  make  your  program  do  the  same  for  the  Hastings  35 
approximation?  Or  do  you  have  to? 


Problem  3 (OPTIONAL) 

Change  your  program  to  use  the  Hastings  Sheet  36  for- 
mula: 


ARCSIN  (X)*ir/2-  (>fT^X)*(l.  57078736-0. 21412453*X+0. 08466649*X*X 
-0 . 03575663*X*X*X+0 . 00864884*X*X*X*X) 


Use  w«3. 14159265 


Problem  4 (OPTIONAL) 

Run  the  program  for  Problem  3 again,  using  Hastings’  37  (on 
cover)  and  the  Double  Precision  feature  of  NBS.  (See  NBS  Primer, 
Section  7,  "EXTENDED  DATA  TiPES".  Use  ir=3 . 141592653589793238  ) 


Problem  5 


In  practice,  the  library  functions  of  NBS  are  used  for 
calculation.  Practice  using  the  ARCSIN (X),  the  ARCCOS (X) , 
and  the  ATAN(X)  functions* in  NBS  to  find  the  angles  A,  B,  C, 
D,  E,  F in  the  figure  below.  (Use  direct  mode  if  you  wish.) 


t 


k* 


Trigonometric  Functions  & 
Tchebychev  Approximations 

(OR — CAN  A HUMBLE  COMPUTER  THAT  ONLY  KNOWS  HOW  TO 
DO  ARITHMETIC  CALCULATE  TRANSCENDENTAL  FUNCTIONS?) 


Problem  1.  Below  you  see  three  representations  of  the  "sine"  func- 
tionT  Which  is  the  ri^it  one? 

□ 1.  (a) 

□ 2.  (b)  . 

□ 3.  (c) 

□ 4.  All  of  the  above. 

□ 5.  None  of  the  above. 


(a)  The  SINE 

function : 


(b)  The  SINE 

function : 


(c)  The  SINE 

functions 


Y 

Y/R 

X 

±/(R2-X2)/R 

A 

SIN(A) 

0.0 

0.000 

1.0 

0.000 

0.0 

0.000 

C.l 

0.100 

0.9 

0.436 

0.1 

0.099 

0.2 

0.200 

0.8 

0.600 

- . 0.2 

0.198 

0.3 

0.300 

0.7 

0.714 

0.3 

0.295 

0.4 

0.400 

0.6 

0.800 

0.4 

0.389 

etc . etc . etc . 


Notes  We  are  using  the 
letters  X,  Y,  A,  and  R 
to  represent  the  same 
quantities  as  in  the  mo- 
dule "Circular  Functions". 
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* 

* TCHEBYCHEV 

Answer  to  Problem  1:  All  of  the  representations  are  correct. 

Here  are  three  NBS  programs  which  you  can  try  to  convince  your- 
self that  you  can  get  three  different- looking  "graphs”  of  the  sine 
function.  The  shape  of  these  graphs  depends  on  which  variable  (X, 

Y,  or  A)  you  decide  to  increment  in  equal  steps. 

After  looking  over  these  programs,  you  may  get  the  feeling  that 
something  has  been  'put  over*  on  you.  Question:  Has  a little  mathe- 
matical "hanky-panky"  crept  into  these  demonstration  programs???? 
Which  one???  (For  the  answer,  read  onlit) 


Demonstration  Program  (a) 


2 


2 


+2 
+ 2 
+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 2 

2+ 

+ 


2 


2 


200 
210 
220 
2 30 
2 40 
2 50 
2 60 
270 
2 80 
290 

295 

296 
2 99 


F0R  Y-0  STEP  *1  UNTIL  Y>*95 
PRINT  TABOO)  :*'+**:  TAB(30+30*Y>  s "2" 
NEXT  Y 

F0R  Y=1  STEP  -*1  UNTIL  Y<*05 
PRINT  TABOO)  9 "+":TABO0+30*Y> 

NEXT  Y 

F0R  Y=0  STEP  -*1  UNTIL  Y<-*95 
PRINT  TAB(30+30*Y):*,2":TABO0 )?•+•• 
NEXT  Y 

F0R  Ys“ 1 STEP  *1  UNTIL  Y>*05 
PRINT  TAB(30+30*Y)  :"2":TAB(30)  :*'+*' 
NEXT  Y 
END 


- + 
+ 

2+ 


# 
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Demonstration  Program  (b) 

+ 1 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

+ . 

i 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

A 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

+ 

i 

1 + 

_ 

I + 

I + 

1 + 

I + 

1 + 

I + 

100  F0R  X*1  STEP  -*1  UNTIL  X<--95 

I + 

110  LET  S*SGRT< 1-X*X> 

+ 

IPO  PRINT  TAB<30):,,+,,:TAB<36+30*S):,,I” 

+ 

130  NEXT  X 

* 

140  F0R  X=-l  STEP  *1  UNTIL  X>-95 

+ 

150  LET  S*-SQRT<1-X*X> 

+ 

160  PRINT  TAB<30+30*S)  J-*M,,JTAB<20) 

I + 

170  NEXT  X 

I + 

■ 

180  END 

I + 

1 + 
1 + 
I + 


14 
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Demonstration  Program  (c) 


300  F0R  A*0  STEP  6.283/40  UNTIL  A>6.283 
310  LET  S=SIN( A) 

320  IF  S>0  THEN  PR.  TAB(  30>  : TAB( 30+30+S)  : "3"  ELSE  PR. 
TAB< 30+30+S) : "3": TAB( 30>  * 

330  NEXT  A 
340  END 


+ 

+ 

+ 3 

+ 3 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 

3 + 


3 


3 


3 


3 


3 + 

3 + 

3+ 


3 


3 


3 


3 


3 


3 


If  you  haven't  spotted  the  hanky-panky  by  now,  check  line  310 
above.  Assuming  computers  can  only  do  arithmetic  (which  is  true), 
an  expression  like  30  + 30  * S in  line  320  makes  sense,  but  it 
isn't  clear  that  arithmetic  is  being  used  in  line  310*  Where  is 
the  arithmetic  being  done?  Answer:  SECRETLY!  The  computer  actually 
goes  into  a "library"  sub-routine  at  line  310,  which  tells  it  what 
arithmetic  operation  must  be  done  on  the  angle  A (expressed  in 
radians)  in  order  to  calculate  SIN  (A) • 


non 


TCHEBYCHEV 


5 


The  Secret  Unveiled 


The  library  functions  in  digital  computers  used  to  evaluate 
complicated  functions  almost  always  employ  expertly  designed  poly- 
nomials as  approximations  to  these  functions.  These  approximations 
can  be  made  quite  accurate,  using  techniques  originated  by  the 
Russian  mathematician  Tchebychev  (also  spelled  Chebyshev) , and  put 
into  useful  form  by  Hastings. ^ 

Here  is  one  of  the  approximations  given  by  Hastings: 

(Hastings  Sheet  14) 

Function:  sin 

Range:  -1  < X < 1 


Approximation : 


* cxx  + c3x* 


where  Cx  = 1.5706268,  C3  = -.6432292,  and  C5  = .0727102. 


Before  using  this  approximation,  it  will  be  useful  to  make  two 
changes.  First,  we  will  change  the  range  of  the  independent  variable 
by  the  substitution: 

Y = (it/2)  * X so  that  when:  -1  < X < 1 

we  will  have:  -ir/2  < Y < ir/2 


Note:  Both  Y and  X are  angles , measured  in  radians. 

. Don't  confuse  these  variables  with  the  Y and  X 
coordinates  used  on  the  cover.  The  variable 
Y in  the  rest  of  this  discussion  is  identical 
to  the  variable  A shown  on  the  cover. 


After  making  this  substitution,  the  Hastings  Sheet  14  approxi- 
mation becomes: 

sin (Y)  = KjY  + K3Y3  + KgY5 

where  Kx  - (2/ir)  * K3  - (23/ir3)  * C3;  K5  = (25/*5)  * Cg. 


la  book  your  library  ought  to  have  is  Hastings,  Cecil,  Approximati ons 
for  Digital  Computers,  Princeton  University  Press,  1955. 


TCHEBYCHEV 
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The  second  change  we  will  make  is  to  rearrange  the  right  side 
of  the  preceding  formula  into  "nested  multiplication"  form: 

sin  (Y)  * (((K5  * Y * Y + K3)  * Y)  * Y + * Y 

The  advantage  of  this  form  is  that  it  only  takes  five  multiplica- 
tions, whereas  the  original  polynomial  takes  nine. 

Here  is  an  NBS  program  which  uses  this  form  to  calculate 
sin(Y),  and  compares  it  with  the  NBS  library  routine  (which  uses  a 
higher  degree  approximation — see  problem  3) . 


PROGRAM  T0  COMPARE  NBS  LIBRARY  SINE  FUNCTION  WITH 
HASTINGS  14. 


RADIANS 

NBS  LIBRARY 

HASTINGS  14 

-I 

-0.84 I 470985 

-0.841534239 

-0.9 

-0.78332691 

-0.783406887 

-0.8 

-0.717356091 

-0.717433028 

-0.7 

-0.644217687 

-0.644277689 

-0.6 

-0.564642473 

-0.564678888 

-0.5 

-0*479425539 

-0.47943851 

-0.4 

-0*389418342 

-0.3894I3I82 

-0.3 

-0.295520207 

-0.295505151 

-0.2 

-0. I 98669331 

-0. 19865316 

-0.  1 

-0.099833417 

-0.099823323 

0 

0 

0 

0.1 

0.099833417 

0.099823323 

0.2 

0.198669331 

0.19865316  . 

0.3 

0.295520207 

0.295505151 

0.4 

0.389418342 

0.389413182 

0.5 

0.479425539 

0.479438509 

0.6 

0*564642473 

0.564678888 

0.7 

0.644217687 

0*644277689 

0.8 

0.717356091 

0.717433028 

0.9 

0.78332691 

0.783406887 

I 

0.841470985 

0.841534239 

>LISTNH 

290  PR.' 

"PROGRAM  TO  COMPARE  NBS  LIBRARY 

SINE  FUNCTION  WITH 

291  PR.' 

"HASTINGS  14.**  PR.  PR. 

300  LET 

P=3. 14 I 59265 

310  LET 

P2-P*P  LET  P3*P*P2 

320  LET 

Cl*( 1 . 5706268*2 )/P 

330  LET 

C3=< -.6432292*8 >/P3 

340  LET 

C5*<. 0727 102*32 )/<P2*P3) 

350  PR. 

DO 

“RADI ANS“*TAB<20)| "NBS  LIBRARY1 

"fTABC 40) S “HASTINGS 

360  FOR  Y*-l  TO  1 STEP  0.1 

365  IF  ABS(Y)<*0001  LET  Y*0 

370  S1=<<<C5*Y*Y+C3)*Y)*Y+C1)*Y 

380  PR.  *'  "SYS  TABC20)  tSIN(Y)  |TAB(40)SS1 

390  NEXT  Y 

400  ENO 


1 


ERIC 


♦ 

y 


f.r 
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Problem  2.  Change  the  FOR  loop  in  the  above  program  to  cover  the 
range  for  Y = 1 to  3 radians#  and  see  what  happens « Also  try 
the  range  for  Y = -1  to  -3.  If  something  ”blows  up,”  can  you 
modify  your  program  so  that  any  value  of  Y can  be  used  in  the 
Hastings'  approximation?  (Hint:  recall  sin(Y)  = sin(ir-Y) 

= sin(2ir+Y)  = sin(3ir-Y)  = etc.) 

. Problem  3.  Here  is  the  Hastings  Sheet  15  approximation  which 

ought  to  come  a lot  closer  to  the  values  produced  by  the  NBS 
library  function.  Modify  your  program  to  try  it  out#  again 
comparing  it  with  the  NBS  library  function. 

sin  ( (it/2)  * X)  ^ Cx  * X + C3  * X3  + C5  * X5  + C?  * X7 

where 

C±  - 1.570794852  C3  - -0.645920978 

C5  = 0.079487663  C?  = -0.004362476 

***  -1  < X < 1 

Problem  4.  (Optional)  If  you  find  the  idea  of  super-precision  of 
interest,  ask  your  teacher  to  see  the  Com- Share  XT RAN  Library 
Manual#  and  explore  the  use  of  double  precision  arithmetic 
and  double  precision  library  functions.  One  of  the  best  hand- 
books to  use  as  a standard  for  "correct"  values  is  "The  Hand- 
book of  Mathematical  Functions#  Applied  Math  Series  No.  55" 
available  from  the  Superintendent  of  Documents#  U.S.  Government 
Printing  Office#  Washington#  D.  c.  20402  ($6.50).  * A paperback 
version  published  by  Dover  Co.  may  be  available  at  book  stores 
for  about  $5.00. 

Problem  5.  (Easy)  Use  the  NBS  library  functions  TAN  and  SQRT  to 
find  the  lengths  of  the  guy-wires  X and  Y shown  below.  (Use 
of  direct  mode  NBS  as  a desk  calculator  is  a good  way  to 
solve  this  problem.)  * / / 


/ 


Phantom  VORTAC 

0 This  module  will  guide  you  in  preparing  the  master  program 
for  an  on-board  flight  computer.  The  computer  takes  informa- 
tion from  a VHF  omni  range  system  (VORTAC) , and  calculates  the 
magnetic  course  and  distance  to  a given  airport  for  the  pilot. 
The  output  of  the  program  is  such  that  the  pilot  can  fly  to- 
ward a "phantom”  VORTAC  located  at  any  airport  he  selects. 

9 A description  of  this  newly  developed  navigational  system 
and  the  mathematics  on .which  it  is  based  are  given  on  page  2. 
9 Pages  3 and  4 outline;  alternate  methods  for  handling  the 
computation^  and  suggest  ways  in  which  previous  programs  you 
have  written  might  be  incorporated  as  sub-routines. 

9 A "real  time"  simulation  of  a flight  using  this  system  is 
suggested  as  an  advanced  level  program. 


V 
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Pilots  flying  over  the  United  States  (and  most  other  countries 
of  the  world)  rely  on  radio  facilities  called  VORTACs  for  navi- 
gational information.  The  basic  information  the  pilot  receives  in 
the  cockpit  is  his  position  relative  to  the  VORTAC,  given  in  polar 
coordinates . 

The  pilot  in  the  illustration  below  would  describe  his  position 
(obtained  from  his  radio  instruments)  as  being. "on  the  100°  radial 
of  the  Allegheny  County  (AGO)  VOR,  50  miles  out". 

It  is  easy  for  this  pilot  to 

note  that  he  can  get  to  AGO  by 

* ** 

turning  right,  and  flying  a course 
of  280° (Why  280°?).  If  he  is 
going  200  miles  per  hour,  it  is  m 
also  easy  for  him  to  estimate  2»o 

that  he  will  arrive  at  AGO  in 
15  minutes  (Is  this  exactly  true? 

— go  back  over  the  module  on  vec- 
tor addition  if  you  are  not  sure.) 

The  catch  to  all  of  this  is 
that  the  location  of  the  VOR  is 
usually  different  from  the  location  of  the  destination  airport. 

This  difficulty  can  be  handled  by  flying  tfro  legs,  the  first  from 
the  present  position  to  the  VOR,  and  the  second  from  the  VOR  to  the 
airport*  This  is  obviously  an  inefficient  route. 


V 


A new  navigational  system 
uses  an  on-board,  special-purpose 
computer  to  tell  the  pilot  what 
course  and  distance  to  fly  in  or- 
der to  go  directly  from  his  present 
position  to  the  airport.  Let's 
first  examine  three  ways  of  analyzing 
the  mathematics  involved  in  such  a 
computation . 


Very  high  frequency  Omni  Range  and  TACAN,  where  TACAN  is  an  older 
military  system  of  distance  measuring  equipment.  Most  civilian 
pilots  call  these  facilities  VOR-DME  stations. 

**The  angular  direction  of  the  intended  flight  path,  measured  clock- 


wise from  N. 


■*9~K. 


VORTAC- 3 


ERIC 


Before  writing  a program  for  such  an  on-board  computer,  it 
would  be  useful  for  a programmer  to  solve  by  hand  typical  problems 

At  this  point,  let's  look  at  three  such 


arising  in  this  situation, 
solutions. 


D=210° 


VORTAC 

E=135° 


Problem  1 . The  aircraft  is 
150  miles  out  on  the  135°  radial 
hnd  the  destination  airport  is 
100  miles  out  on  the  210°  radial. 
That  is,  in  the  diagram.  Si  = 150, 

E = 135°,  S3  = 100,  D = 210°. 

Let  A2  be  the  angle  determined 
by  Si  and  S3.  Hence,  A2  = 

210°  - 135°  = 75°.  By  the  Law  of 
Cosines, 


S2  = /1002  + 1502  - 2 


100 


150  * cos (75° ) = 157.3  miles. 

How  by  the  Law  of  Sines,  sin  A3  =*  (100  a .966)/157.3  « .614  and 
cos  A3  = /l  -,6142=  .789.  From  this  we  can  get  A3  = arctan(.614/.789) 
= 38°.  (Why  wasn't  A3  computed  directly  from  sin  A3?)  Therefore 
X ss  135°  + 180°  - 38°  = 277°.  Output  to  the  pilot  is  157.3  miles, 
277°— is  this  answer  reasonable? 

Problem  2.  The  airplane  is  250  miles  out  on  the  45°  radial 
and  the  destination  airport  is  100  miles  out  on  the  280°  radial, 
in  terms  of  the  diagram.  Si  - 250,  E = 45°,  S3 1 * 100,  D « 280°. 

Converting  the  positions  of  the  aircraft  and  destination  to 
rectangular  coordinates  with  the  VORTAC  as  origin  we  have 

sin  45°  = 176.8  miles, 
cos  45°  = 176.8  miles, 
sin  280°*  -98.48  miles, 
cos  280°*  17.36  miles. 


x - 250 

P 

y = 250 
P 

= 100 


100 


A = arctan 


176.8  - 17.36 

176.8  - (-98.48) 
arctan  (.5792)  = 30.1° 


VORTAC 


This  is  another  way  of  saying  that  computers  do  not  remove  the  re- 
sponsibility of  analyzing  a problem  before  solving  it — quite  the 
contrary;  they  demand  more  thought  than  ever.  This  may  be,  in  fact, 
one  of  the  most  important  contributions  computing  systems  can  make 
to  learning. 


22 
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Therefore  X = 270°  - 30 . 1°  = 239 .9°  . By  the  Pythagorean  Theorem 
S2  = / (176 . 8~-  17. 36)2  + (176.8  - (-98. 48))2  = 318.1.  Hence  the 
output  to  the  pilot  is:  318.1  miles,  239.9°. 

Problem  3.  The  aircraft  is  150  miles 
out  on  the  120°  radial  and  the  destination 
airport  is  400  miles  out  the  150°  ra- 
dial.** In  the  diagram,  £ - 120°,  Si  = 150 
miles,  D - 150°,  S3  - 400  miles.  Using 
the  same  reasoning  as  in  Problem  1 yields: 

A2  = 30° 

S2  = /4002  + 1502-  2 • 150  • 400  • cos  30° 

= 280.3  miles 

sin  A3  = 400  . (.5236/280.3)  « .7135 
cos  A3  = /I  - .7135 2 = .7007 
A3  = arctan  (.7135/. 7007)  = 45.52° 

X * 120°  + 180°  - 45.52°  = 254.48° 

Output  to  the  pilot  is:  280.3  miles, 

254.48° — does  this  answer  look  reason- 
able? What  went  wrong? 

Assignment.  Write  a program  for  an 

on-board  navigational  computer  which  will 

accept  as  input  the  radial  and  distance  from  a VORTAC  of  both  an 

aircraft  and  a destination,  and  which  will  compute  a course  and 

distance  to  the  destination.  The  trigonometric  subroutines  in  the 

computer  require  arguments  in  radians,  but  pilots  think  in  terms 

of  degrees,  so  it  will  be  necessary  for  you  to  convert  degrees  to 

radians  and  back  (see  the  module  on  converting  to  radians) . Other 

possibly  useful  topics  are  inverse  trigonometric  functions  , the 

* 

Laws  of  Sines  and  Cosines  , and  transformation  of  polar  to  rectan- 

„ * _ - 
gular  coordinates  . 


* 

See  modules  you  have  worked  previously. 

** 

By  now  you  should  have  noticed  that  the  word  radial  is  used  to 
designate  the  angular  position  of  a line  segment  that  starts  at  the 
VORTAC. 


N 

1 
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Optional  Section.  During  most  flights,  the  aircraft  moves 
through  the  air  and  the  air  moves  as  well.  Write  an  addition  to 
your  program  which  will: 

1.  Accept  as  additional  input,  specified  by  the  'pilot', 

(a)  the  aircraft's  speed,  (b)  the  speed  of  the  wind,  (c)  the  direc- 
tion of  the  wind,  (d)  an  elapsed  time,  t,  and  (e)  a heading* for  the 
aircraft.  Any  heading  should  be  acceptable:  the  'pilot'  should  be 
able  to  fly  wherever  he  likes. 

2 . Compute  the  new  position  of  the  aircraft  on  the  basis 
of  the  above  information. 

3.  Repeat  steps  1 and  2 as  often  as  desired  or  necessary  to 
reach  the  destination  airport. 

You  can  look  at  the  input/output  sequence  of  one  program  that 
accomplishes  this  and  has  been  written  and  stored. 

In  order  to  gain  access  to  this  program,  please  logon  (direc- 
tions are  given  in  an  introductory  module)  using  the  User  Identifi- 
cation Code  given  to  you  by  your  teacher: 


t; 

s 


In  order  to  interact  with  this  stored  program  type  in  the  statements: 


Prom  this  point  on  the  operation  of  the  program  should  be  self- 
explanatory. 


* 


Heading  is  defined  as  the  angular  direction  of  the  longitudinal 
axis  of  the  aircraft  with  respect  to  North.  In  the  picture  below, 
the  pilot  is  flying  a course  of  90°,  but  his  heading  is  65°.  Why? 
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GROSS  FLOW  CHART  FOR  SUGGESTED  ASSIGNMENTS 


Aircraft 

Position 

Airport 

Position 


